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We show that doping-induced charge fluctuations in strongly correlated Hubbard electron sys- 
tems near the |-fllled, insulating limit cause overscreening of the electron-electron Coulomb repul- 
sion. The resulting attractive screened interaction potential supports dj.2_y2-pairing with a strongly 
peaked, doping dependent pairing strength at lower doping, followed by s-wave pairing at larger 
doping levels. 
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In the cuprate high- Tc superconductors, the strong, 
on-site Hubbard- J7 Coulomb repulsion is believed to pre- 
vent conventional on-site Cooper pair formation. This 
long-standing theoretical dictum seems to be now sup- 
ported by experimental evidence, obtained in several 
cuprate materials, for a non-s-wave pairing state, appar- 
ently of d2;2_y2 -symmetry 



We start from an extended Hubbard Hamiltonian 
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Ht , 
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terial, the {Nd, Ce)2CM04-system, appears to exhibit a 
fully gapped excitation spectrum, suggesting an s-wave 
state within the framework of more conventional pair- 
ing theories A d-wave pairing state implies a spa- 
tially extended Cooper pair wavefunction involving elec- 
trons paired at 1st neighbor or larger lattice distances. 
While impervious to the on-site U repulsion, even non- 
s-wave extended pairing states can be suppressed by the 
extended (1st, 2nd, ... neighbor) part of the electron- 
electron Coulomb repulsion. For any proposed micro- 
scopic model of the cuprates, it is therefore of crucial im- 
portance to demonstrate that, firstly, the model can sup- 
port at least two qualitatively different pairing states of, 
possibly, different symmetry Q and, secondly, that such 
pairing states are robust against on-site and extended 
Coulomb repulsions. 

In the present paper, we use a combination of dia- 
grammatic and quantum Monte Carlo (QMC) techniques 
to show that the charge fluctuations in a quasi two- 
dimensional (2D) extended Hubbard model give rise to 
a screening effect which not only reduces the magnitude, 
but in fact reverses the sign of the extended part of the 
three-dimensional (3D) screened Coulomb potential Vs- 
At larger doping, this overscreening becomes so strong 
that even the on-site part of Vs changes sign and turns 
attractive. The overscreening is intrinsically a large- C/ 
effect of the Hubbard system at finite doping near band- 
fiUing i. Over a doping range of about 10 — 20%, \Vs\ is 
much weaker than the bare U and thus offers the prospect 
of developing controlled perturbative expansions in terms 
of V5, rather than U. Taken as an effective pairing poten- 
tial, Vs supports dr^2_y2 pairing, with a pairing strength 
maximum in the ~ 10 — 20% doping range, and, in close 
proximity, s-wave pairing at larger doping. 



Yet, at least one ma- with c] ^ creating an electron of spin a =1, J. at Cu-site rj 



in a 3D crystal of stacked Cu02 layers, nj = cjcrCjo- 
and rj£ — rj —r(. Ht includes only an in-plane 1st neigh- 
bor hybridization t and the chemical potential /i. The 
3D Coulomb potential 

Vir) = USr.O + -^{l - Sr.o) = USr.O + K(r) (2) 

includes the on-site (r — 0) repulsion U and an extended 
l/|r|-part, T4, with a dielectric constant es to account 
for screening by the insulating background not explicitly 
included in H, that is " non-Hubbard" electrons in lower 
filled bands and, possibly, phonon degrees of freedom. 

The basic idea of our approach is to treat some short- 
range portion of V{r), denoted by Vo{r), exactly by QMC 
methods. The remaining weaker, but long-range part of 
V, denoted by Ve{r), is then handled perturbatively. By 
including only short-range in-plane terms in Vo, the QMC 
simulation can be restricted to a single 2D layer for which 
we obtain, by QMC, the density correlation function 
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-Po{q,luj)[l-Vo{q)Po{q,lu:)]-^ 



(3) 



at wavevectors q, Matsubara frequencies ioj and temper- 



ature T = 1/(3 for lattice size N with Anj 



Here, (...)o and ...(t)o denote, respectively, thermal aver- 
aging and imaginary-time evolution with respect to the 
QMC Hamiltonian Ho = Hy^ +Ht. Pq is the correspond- 
ing irreducible polarization insertion Q . 

The exact screened potential Vs, irreducible polariza- 
tion insertion P and density correlation function x of the 
full Hamiltonian (|]) are related by ||^ 
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Vs{q,iu;)= [l-V{q)P{q,iu;)]-'V{q) 
= V{q)~Viq)x{q,tto)Viq) , 



(4) 



where V{q) denotes the lattice Fourier sum over V{r). 
Our essential approximation is to replace the exact P of 
the full Hamiltonian H hy Pq, extracted from the QMC 
results for Xo via Eq.||. All renormalizations of P due to 
Vo are thus included exactly, to all orders in Vo- Renor- 
malizations due to the weaker long-range part Vi are ne- 
glected in P, but approximately included in Vg via Eq. ^. 
Note here that Eqs. |3[ ^ for V5, P and Pg are based on 
the diagrammatic expansion in the charge representation 
where the /7-term is written as ^ nj, rather than the 
more familiar spin representation UY^jUjiriji. While 
both are equivalent when summed exactly to all orders, 
the former, as we will show, has some crucial advantages 
for approximate diagram resummations. 

As a simple cuprate, we consider La2-xSrxCuOji, with 
a body-centered tetragonal {hot) model crystal structure, 
in-plane lattice constant a — 3.80A, inter-layer spacing 
d = 6.621 %,t = 0.35ey and t/ = 8t §. Using standard 
determinant QMC methods, Xo is simulated with up to 
2 X 10^ MC sweeps and typically < 0.5% statistical error 
on 6 X 6, 8x8 and 10 x 10 2D lattices with periodic 
boundary conditions, taking Vo{r) = USrfl, at /? = l/T 
up to pt — 10, with imaginary-time step At = < 
0.0625^^^ where L is the Trotter number. 

To estimate es from the measured long-wavelength ex- 
ternal field dielectric function eea;(w)|(;— »o in the undoped, 
insulating La2CuOA &0| we calculate Q 
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at the smallest q- vectors with non-zero in-plane compo- 
nent available on our finite model lattices and iuo = Q, 
with Vc denoting the 3D unit cell volume, x = ~P[^ ~ 
V{q)P]-^ from (|) and P ^ Po- In the undoped Hub- 
bard system at temperatures T > O.lt, we find that eex 
from (^) varies roughly linearly with e^. The excess 
Ae = Cea; — e_B — 0.75 — 0.85, that is, the Hubbard elec- 
trons' contribution to the dielectric screening, is approxi- 
mately independent of eg for es > 3. From the measured 
dielectric function e^x of undoped La2CuO^ in the static 
limit ^ 0, eo = texi^) — 30 ^, we thus estimate 
eb = Cex — Ae = 29. However, this includes a large, 
in fact, dominant phonon contribution pl| . The purely 
electronic dielectric screening is observed at frequencies 
Woo ~ 0.5 — leV, well above the phonon spectrum flph < 
O.leV, but still well below the electronic Mott-Hubbard 
charge gap Amh 1.5 - 2eV, where Coo = £ex{uJoo) = 5 
1^. Hence we get €b =4.2 for the purely electronic back- 
ground screening. /,Prom the estimated values e_B =4.2 
(without phonons) or even es = 29 (including phonons), 
one obtains a quite substantial 1st neighbor repulsion 
strength Vi = e^/(eBa) = 0.90eV = 2.6t in the former 



and Vi ^ 0.13eV ^ 0.37^ in the latter case. Thus, Ve(r) 
could severely suppress spatially extended pairing poten- 
tials (l|. 

In Fig.|l|(a) and (b), we explore how screening affects 
the 3D Coulomb potential Vs, calculated from (^) at 
ioj — 0, Fourier transformed back to r-space and plot- 
ted vs. doping concentration x = 1 — {nj) for the on-site 
(r = 0) and in-plane 1st neighbor r- vector. The surpris- 
ing result in Fig. ^(b) is that a small amount of dop- 
ing, of order 5%, will not only suppress the extended 
l/|r|-repulsion for r ^ 0, but will in fact cause a sign 
change in the 1st neighbor and similarly (not shown) 
in the 2nd and 3rd neighbor screened potential. Thus 
Vs{r) becomes attractive at short-range distances. The 
attraction strength at r ^ shows a pronounced doping 
dependence, reaches a maximum at x ^ 10 — 14% and 
becomes repulsive again at x~23 — 28%. Even more sur- 
prising, as shown in Fig. |l](a), is the doping dependence 
of the on-site (r = 0) potential which is also rapidly sup- 
pressed with increasing x and becomes strongly attrac- 
tive at larger doping, near x= 15%. Minus sign problems 
at finite doping unfortunately limit our simulations to 
T>0.33t. However, at least in that temperature regime, 
we find IV5I to be increasing with decreasing T. This 
suggests that the overscreening becomes stronger than 
shown in Fig. |l| at lower T. 

The presence of a strong Hubbard-J7 and finite doping 
density a; > are crucial for the overscreening. If one 
replaces the V"o-renormalized Po in (^ by, say, the non- 
interacting (" RPA" ) polarization bubble Prpa , one also 
obtains a suppression of Vs{r). However, both the on-site 
(r = 0) and the short-range extended part (r 7^ 0) of Vs 
remain repulsive in RPA p^ . Likewise, in the undoped 
large-C/ system, Vs{r) is only reduced relative to V{r), by 
a roughly r-independent factor, comparable to the ratio 
<^ex/^B, for r 7^ 0. Thus, Vs retains a l/|r|-dependence 
and remains repulsive. This is expected for the screening 
of a Coulomb potential in an insulator and confirms the 
insulating character of the i-filled Hubbard system. 

To understand the central role of U and finite doping x 
in the overscreening effect, consider the pure 2D Hubbard 
model where V{r) = Vo{r) = USrfl, our approximation 
P = Po becomes exact and, from Eq. 0, 



Vs{r,iu!) = USrfi - U^Xo{r,iLo) 



(6) 



with Xo(^:*'^) denoting the Fourier transform of 
Xo{q,ii^)- Clearly, the on-site potential Vs{r = 0,iuj) at 
finite doping must become attractive for J7 — > 00, since 
XoiT — Ojiuj) is always positive and approaches a non- 
zero independent limit [of, at least, 0{x/t), by a sim- 
ple U = 00 scaling argument] for U ^ t. By contrast, 
at i -filling, all charge fluctuations are suppressed, with 
Xo{r, = 0) ~ 0{t^/U^), and Vsir = 0,iuj = 0) ^ U is 
repulsive for U^t. For near (1st, 2nd, ...) neighbor r's, 
Xo [f^ iuj = 0) is negative at small U and becomes posi- 
tive only at finite doping and only when U exceeds some 
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FIG. 1. Screened Coulomb potential Vs(r) = Vs(r', ioj = 0) 
at (a) on site (r = 0) and (b) at in-plane 1st neighbor lat- 
tice vectors r and (c) Eliashberg A-parameters for in-plane 
1st neighbor d^-2_yi and on-site/lst neighbor s-wave pairing, 
all plotted vs. hole doping concentration x = 1 — {rij} at 
l3t = t/T = 3.0 and Art = 0.0375 for es = 4 and 30. In 
(a) and (b), results are for La2Cu04, with bet lattice sizes 
6x6x6, 8x8x8 and 10 x 10 x 10, with estimated statisti- 
cal uncertainties in Vs of less than O.OleV. In (c), results are 
based on 8 x 8 x 8 data for Vs- Xd is scaled x4 for display. 



doping dependent threshold of order several t. The near 
neighbor attraction in Vs therefore also requires large U 
and finite doping. 

The results in Fig. |^ are almost independent of Ve for 

>4. (The = 30 data are within 1 — 2% of the pure 
Hubbard (es — > 00) results.) By using spatially trun- 
cated versions of Ve, different layered geometries, and 
more realistic dielectric background models [|lO| , we have 
also verified that neither 3D Coulomb interlayer terms, 
nor the long-range l/|rj-tail in Ve, nor our particular 
choice of crystal structure, nor local field effects will sub- 
stantially affect the overscreening. To test our approxi- 
mation, P = Po, we have also carried out QMC simula- 
tions with both U- and a 1st neighbor Vg-term included 
in Vo- Preliminary results suggest that setting P = Po 
in (Q) reproduces the main T4-effect on x, but tends to 
underestimate the attraction in Vs- Thus, the essential 
features of Vs are very robust against extended Coulomb 
terms. This result can be rationalized by expanding (Q) 
to leading order in Ve- The corrections are smaller than 
the C/^Xo-term in (|6|) by a factor of order ViXo for U^t 
and finite x^t/U- By simple U—^oo scaling arguments 
this is of order Vix/t, i.e. small ii Vi t/x. At suf- 
ficiently large Vi, i.e. low es, our basic approach does 
break down due to charge density wave instabilities, sig- 
naled by l/x(9: iuj = 0) —>■ 0. For the present parameter 
set and doping range, this happens only for es < 1.8. 

A crucial advantage of our diagrammatic expansion in 
the charge representation is the large reduction of the 
overall strength of Vs in the 10 — 20% doping range. This 
suggests the possibility of carrying out controlled, self- 
consistent weak-coupling expansions in which the fully 
screened Vs, rather than the bare V or U, serves as the 
small parameter. As a first step in that direction, we have 
explored possible T^s-induced or -enhanced superconduct- 
ing pairing instabilities, using the standard Eliashberg- 
McMillan (EM) approach |Q. A convenient measure of 
the pairing strength of Vs are the dimensionless EM A- 
parameters, defined in terms of the Fermi surface (FS) 
"expectation values" of Vs{k — k',iu! = 0) for relevant 
Cooper pair trial wavefunctions ri{k) in electron momen- 
tum (fc-) space, as described in detail in Refs. fl^ . 

In Fig. 1(c), we show the EM parameters As for on- 
site s-wave (and, identically, for in-plane 1st neighbor 
s-wave), and A^, for in-plane 1st neighbor dx2_y2 pair- 
ing, with respective pair wavefunctions T]s{k) = 1 and 
Vdik) = cos{akx) — cos{aky). For the required FS in- 
tegrals, we have interpolated our 3D Vs{q,iuj) from the 
finite-lattice 8x8x8 onto a 200 x 200 x 200 g-mesh. At 
low doping, the dominant attractive (A > 0) channel is 
dx2_y2 with Arf reaching a maximum of ^0.15 — 0.17 near 
10 — 14%. As is repulsive at low doping, but becomes 
strongly attractive at larger doping x > 15%. Thus, as 
expected on symmetry grounds. As and A^ reflect the 
doping dependence of the on-site and 1st neighbor at- 
traction Vs shown in Fig. |^(a) and (b), respectively. The 
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A-values for near-neighbor pair wavefunctions of other 
symmetries {p, d^y, g) are smah compared to As and A^. 

The spectral weight of xil^^^) extends up to values 
^^x~8 - lot In the EM analysis Hi, this "boson" 
energy scale, together with A, determines the supercon- 
ducting Tc, roughly as Tc ~ fi^ exp(— 1/A). Because of 
the large fi^^-scale, it may be possible to achieve high T^s 
even at moderate coupling values A< 1. 

It is quite possible that near- neighbor, c?- wave- 
attractive charge correlations x(r 7^ 0, iw) > are very 
closely related to short-range antiferromagnetic spin cor- 
relations in Hubbard systems near ^-filling. The charge 
fluctuation picture developed here may thus provide a 
description of the physics in near-i-fiUed Hubbard sys- 
tems which is complementary to that of a spin fluctua- 
tion based approach |^ . The overscreening of the on-site 
potential, and hence the possibility of s-wave pairing in 
the Hubbard model, is one aspect of this problem which 
is obvious in the former, but difficult to capture in the 
latter approach. 

In summary, we have studied the effect of screening 
on the electron-electron interaction potential in a quasi- 
2D Hubbard model for CM02-layers, coupled by an ex- 
tended l/\r\ 3D Coulomb repulsion. While contribut- 
ing only a minor portion of the total insulating dielectric 
constant at ^-filling, the Hubbard electron system, when 
doped away from ^ -filling, exhibits a strong overscreen- 
ing effect which causes the extended part of the screened 
potential to change sign and become attractive, at 1st 
and further neighbor distances. This screened potential 
gives rise to a pairing attraction in the d^^-yi channel 
which, as a function of doping, exhibits a maximum near 
'--^ 10 — 15%, reminiscent of the doping dependence of the 
superconducting Tc in the cuprates. At larger doping, 
even the on-site part of the screened potential becomes 
attractive and gives rise to an s-wave pairing attraction 
which increases strongly with doping, suggesting the pos- 
sibility of a doping- induced transisition or cross-over from 
d- to s-wave pairing. The overscreening effect is robust 
against 3D extended Coulomb repulsions, independent 
of the 3D crystal structure, and represents intrinsically a 
charge fluctuation aspect of the Hubbard electron system 
at large U and finite doping density near ^-filling. 
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